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ABSTRACT 

Considering a general mass matrix for quasi-degenerate neutrinos and treating the ex- 
perimentally known oscillation parameters as inputs, we study the correlation between the 
degenerate mass scale (m) and the mixing angle ^13. We find that, corresponding to different 
values of m, there exist upper bounds on 6'i3, so that a precise determination of the latter 
in future may put upper limit on the former, and vice versa. One can also find a possible 
correlation between m and lower bound of ^13, depending on the relative strength of the 
unperturbed degenerate mass matrix and the perturbation. The possible constraints on the 
parameters of few models of quasi-degenerate neutrinos are briefiy discussed. 
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The currently available experimental data on neutrino oscillation indicate that the mass 
squared differences of neutrinos are quite small. Although the mass-squared differences and 
mixing angles (except 613) are known from neutrino oscillation experiments, the absolute 
values of the neutrino masses still remain unknown. Depending on the scale, the pattern 
of neutrino masses may be hierarchical or quasi-degenerate in nature. In particular, if one 
assume that relic neutrinos constitute the hot dark matter of the universe, then the scale of 
neutrino masses are expected to be somewhat higher than their differences, and neutrinos are 
expected to be quasi-degenerate. Considering the recent results ^Lj from WMAP experiment, 
there might be some improvements in the upper bound on the neutrino mass. However, 
keeping in mind all possible uncertainties in the cosmological bounds, an upper limit of 2 
eV on the sum of the masses of three neutrino flavours can be said to exist [2j. 

Quasi-degeneracy of neutrinos has also been suggested in a number of theoretical models 
proposed in the literature [3]. For example, considering neutrino masses as degenerate at 
some seesaw scale, various authors [U El El El El El have shown that large mixing angles 
for solar as well as atmospheric neutrinos can be obtained after extrapolating masses and 
mixing to the weak scale using renormalization group equations. 

Among the neutrino oscillation parameters, the values of two mass-squared differences 
and two mixing angles, viz. 612 and 623, are already known to a reasonable degree from 
the solar llO\ and atmospheric [H] neutrino data. However, quantities still unknown are the 
mixing angle ^13, mass scale m of the neutrinos and the CP violating phase 6 in the neutrino 
mixing matrix. It should be noted that the CHOOZ-Palo- Verde experiments fi2[ [13] have 
put an upper bound on 6*13 (< 12°). Furthermore, the absence of neutrino-less double beta 
decay implies an upper bound on m [l4j, which is compatible with the limit coming from 
the hot dark matter content of the universe [Tj. 

In this note, we would like to illustrate the interplay among the unknown quantities, 
namely, m, elements of the perturbation matrix and 613, assuming that neutrinos are of 
Majorana nature and the masses are quasi-degenerate. We have taken into account the 
solar and the atmospheric data as well as neutrino less double beta decay constraint in our 
analysis, which has not been considered in earlier analyses. As we shall see, the maximum 
allowed value of ^13 and the cosmological upper limit on the overall neutrino mass scale 
actually restrict the parameter space of the perturbation matrix lifting the degeneracy. This 
can enable one to constrain various theoretical models of the neutrino mass matrix, which in 
turn determine the way degeneracy is lifted. We will also try to find any possible correlation 
between the degenerate mass scale m with minimum value of ^13. 

Without going to a specific model, we shall first consider the most general form of de- 
generate mass matrix in the weak interaction basis, which allows mixing among different 
flavours of neutrinos due to different intrinsic CP-properties. This degeneracy is lifted, for 
example, by the breaking of some symmetry at the seesaw scale [15j. Thus we add a small 
perturbation matrix to the original degenerate mass matrix in a model-independent way. 
We use as our inputs the known oscillation parameters, namely: 

• I Am^gl ~ 2.12+^:°5' X 10-3 eV^, 623 ~ 45.0°+^°3p° (from the atmospheric z/^ deficit [llj). 

• Amf^ ^ 7.9^^:^ X 10-5 eV2, 612 ~ 33.21°lti° (from the solar z/^ deficit [lO]). 

where Amfj= m?^ — rnf. 

Thereafter, the application of degenerate perturbation theory, in conjunction with the 
constraints arising on the perturbation matrix after using the above experimental values. 
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allows us to obtain the correlation sought after. Specific choices of models for generating 
neutrino mass are likely to restrict further the general form of the perturbation matrix. We 
study the correlation in the context of a few models having family symmetries of abelian 
and non-abelian nature, and check whether these in turn can constrain some parameters of 
the models even further. 

With the charged lepton matrix taken as diagonal, real and positive, the Majorana mass 
terms in the flavour basis can be expressed as 

ll^mass = -{viSC'^^apVLp + h.C (1) 

where Mq^ is a 3 x 3 complex symmetric mass matrix and is the weak eigenstate basis of 
neutrinos corresponding to three generations. The neutrino flavour states |z/q,), a = e,fi,T, 
in the weak basis are related to the neutrino mass eigenstates i = 1,2,3, with masses 
rrii : 

Wa) = ^ UaiWi) ■ (2) 
i 

where [/ is a 3 x 3 unitary matrix. In general, the mass matrix M can be diagonalised by a 
transformation of the form 

mi 

U'' MU = Mdiag = I m2 I (3) 

mg 

where Mdiag is the diagonal mass matrix. If one considers the three masses to be degenerate, 
then U can be rotated away for Dirac neutrinos. This cannot, however, be done in the case of 
Majorana neutrinos if the CP property of one of the fields is different from those of the other 
two. It was shown by Branco et al [16] that, in such a case, the fact that the CP-eigenvalue 
of a Majorana state can be +i or —i implies a still non-trivial form of the the diagonalising 
matrix U which contains two mixing angles and one phase. Implications of such (or similar) 
scenarios have also been explored, for example, in reference 

For degenerate Majorana neutrinos remembering that the mixing angles 612 and 6*23 are 
required to be large by observation, we shall consider the intrinsic CP-eigenvalue associated 
with z/2 to be opposite to that of the other two neutrinos. Without losing any generality we 
can thus write 



U= \ 43 ^'12 -c'i2 (4) 
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where c-^ = cos6'^, s-^ = sin6'^ and = 7r/2 corresponds to different intrinsic CP property 
of V2 with respect to other two neutrinos and a corresponds to CP violating phase in general. 
We shall choose it to be zero for the CP conserving case. Following ref. the degenerate 
mass matrix M in the flavour basis can be written as 

M = m U*U^ (5) 

where degenerate mass scale m = nii = 1112 = m^. 

However, in spite of such a mixing, neutrinos cannot oscillate, since the quantity governing 
oscillations are the mass-squared differences. In order to lift the degeneracy, we assume that 
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degenerate mass scale m is somewhat higher than the required mass-squared differences as 
indicated by the oscillation experiments. There are studies on the mechanism for lifting such 
degeneracy. One can assume degeneracy at a high scale and envision that it is lifted through 
running, as some symmetry which holds at high scale is broken. Alternatively, one can 
remain confined to the electroweak scale itself and consider mass splitting effects there. In 
either case, the mechanism consists in a perturbation to the mass matrix, in a basis where it 
is diagonal and degenerate. This perturbation matrix is real and symmetric in several studies 
in recent past. We follow the same practice here. Moreover, most of these investigations do 
not throw much light on the numerical values of the elements of the perturbation matrix. 
This is the point we wish to address, namely, how the limits on ^13 as well as the neutrino 
mass scale restrict the values of the perturbation matrix elements. We believe that , in spite 
of a six-fold multiplicity of the elements, such constraints can be useful in shortlisting viable 
models. 

We thus consider a small perturbation to the degenerate mass matrix M^iag, parameter- 
ized as 




Q^el a g 1 \ (6) 



where is of the order of 10~^ eV^ and other parameters are < 1. In general, the smallness 
of present in the standard parametrisation of neutrino mixing matrix (as mentioned in 
equation (15) below) hints at the CP violating effect being small in the neutrino sector. So 
keeping CP violating phase only in M we have neglected it in the small perturbation matrix 
Q making it real. Q is written in a basis in which the degenerate mass matrix is diagonal. 
It is also obvious that in the flavour basis Q becomes 

Qo^U*Q (7) 

where U is of the form indicated in (4) . Unlike Q in the flavor basis Qo is complex. We next 
employ the methodology of degenerate perturbation theory, with the diagonal and degenerate 

mass matrix determining the unperturbed basis. Besides, we absorb the phase associated 
with /9 = 7r/2 in the Majorana neutrino field, thus enabling it to disappear from the matrix 
U and we define this as Ui. Considering M in (5) and using Ui as the diagonalising matrix 
from eq. (3) we get the diagonal elements of M^iag as , 

777.2 = — ?77, mi = 7773 = 777. (8) 

Now, only rrii and 777.3 are degenerate. If we wish to parameterize the lifting of degeneracy 
in such a manner that the mass 777i remains unchanged, then one may set the following 
condition on the parameters in Q: 

= ef (9) 

Using first order degenerate perturbation theory, we obtain the following mass eigenvalues 
for Mdiag + Q, lifting the degeneracy to first order in e: 

777i = 777, 7772 = " ?77. -|- eg, 7773 = 777 -|- e(e -|- /) (10) 
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and the diagonalising matrix as 



U' 



( cosip —ar sinip 

r {a cos ijj — sin ijj) 1 r {a smijj + cos ip) | Af (11) 
^ — sin ip —r cos ip 



where tan 2?/' = 2b/ {f - e); r = J\f = diag{N{\ N^^) where A^i, A^2,A^3 are the 
proper normahsation constants for the eigenvectors. These can be easily expressed in terms 
of the parameters of the theory. 



A/i = [l + r^(a cos ^/^ — sin ^z^)^] ^ 
As = [l + r^(a sin + cos '?/')^] ^ 

(12) 

To go to the physical mass eigenstate basis, however, we keep j3 = tt/2 in U and now 
the eigenvalues are same as in equation [8] except there will be overall change in sign in m2 
making it positive. So replacing Ui by U we write 

Uo = UU' (13) 

which will diagonalise the mass matrix M + Qq. Although there was no mixing of and 
I'r in U but Uo has that mixing. One may note that, since the solar neutrino data indicates 
that 1712 is heavier than mi and that there is no sign ambiguity in Amfg, the parameter g 
should be negative in our case. To relate the different parameters with the experimental 
data we shall write Uq in the standard form of a 3 x 3 unitary mixing matrix: 

C12C13 ■S12C13 Si3e 

Uq = I — C23S12 — S23Si3Ci2e*'^ C23C12 — 5235135126**^ S23C13 | , (14) 
S23S12 - C23Si3Ci2e*'' -S23C12 - C23Si3Si2e*'' C23C13 

where Cij = cos 9ij, Sij = sin 9ij and 6 is the CP violating phase. 

The parameters of Q in equation [6] can be constrained from the mass squared differences 
obtained from solar and atmospheric neutrino data. Using equation [10] up to the order of e, 
one can write 

2meg = Ami2, 2me(e + /) = Am23 (15) 
From equation [9] it follows that e and / cannot be of opposite sign which implies 

Qii =ee< ^m\j2m (16) 

It follows from the above that 

x(= em^) < Am23|/'r (17) 
We consider ^c? < [Am^gl so that solar neutrino data is satisfied. 
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We further simplify the analysis by setting 5 = in equation [H] and a = in equation 
m which correspond to CP conserving case. Comparing 11, 12 and 13 elements of Uq in eqs. 
and [m we obtain 



sin 6* 



13 



cos sin ip + sin {o, sin ip + cos ip) 



2m 



where 



tan 6\ 



12 



aci2e + 2msi2 cos ip 



2mci2 — easi2 cosip + esu sinip 
Equations [18] and [19] can be re-expressed as in term of x (= em^ 



(19) 



sin 9 



13 



and 



cos6'^ 



12 



2rx 



+ r sin ^1 



12 



'1- 



2rx 



+ a* 



2rx , 



(20) 



tan6'^ 



A[2 



tan^^ 



12 



2ra) 
Amis 



+ ar 



12 



r I a, /l — 



_2nr 



2ra: 



(21) 



Considering 612 from the solar neutrino experimental data as discussed in the introduc- 
tion, one can find the relation between the two unknowns m and ^13 in terms of the elements 
of Q. It may be noted here that for both normal and inverted hierarchy equations [20] and 
Un remain unchanged. 

Now one can see how ^13 and m are related. It is evident that for a given value of r, ^13 is 
a function of x only. A careful inspection of equation [TH tells the monotonically increasing 



> 




Figure 1: Allowed region for m and e depending on ^13. r = = 0.1 has been used in this plot 
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dependence of ^13 on x. Thus maximum allowed value of ^13 corresponds to the limiting 
value of X as defined in equation [T71 The contours of sin 6*13 thus correspond to specific 
values of X. In other words, they corresponds to contours of specific values of x in the e — m 
plane. 

In figure [H we vary e in the range from to its upper bound defined by equation [T71 
Corresponding to different values of 613 we have shown the variation of m as obtained 
numerically. The region marked with 613 > 12 is thus disallowed from the CHOOZ result. 
We have argued earlier how the normal and inverted hierarchical ordering of neutrino masses 
would produce the same constant 613 contours in the m — e plane. Depending on the possible 
values of e, one can find an upper bound on m in both normal and inverted hierarchical case 
from this plot. If, in future, better bound on 6*13 is obtained then it is possible to improve 
upper bound on the degenerate mass scale, m, compared to its present cosmological bound 



1-8 I 1 1 1 1 1 r 




I I I I I I I I I 
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Figure 2: Minimum value for ^13 vs. the degenerate mass scale m for two different values of 
perturbation parameter e. 



In figure [2], we have plotted the minimum values for 613 as a function of the degenerate 
mass scale m for two different values for the perturbation parameter e. The range over which 
we vary m is fixed by two considerations. The upper limit should be less than 1 eV from 
cosmological considerations [T] . Furthermore, we would like to generate the tiny mass square 
differences by perturbing the degenerate matrix (with eigenvalues m). Thus the value of m 
should be bigger than the differences (which is generated by perturbation) in the masses. 
Thus, keeping in mind the atmospheric mass (square) difference, we take the lower limit for 
m to be 0.1 eV. 

It is interesting to note that that there is lower bound on 6*13 depending on the value of m. 
One can see that the minimum of 613 decreases with increasing value of m. Increasing m and 
keeping the value of e fixed implies, the relative strength of the perturbation with respect to 
the unperturbed degenerate mass matrix, keeps on decreasing. This in turn is taming down 
the value of ^13. In the limit when perturbation vanishes, masses are exactly the same and 
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^13 is exactly zero. However, to create qusi-degenerate masses there is perturbation and as 
such there is some lower bound on ^13 which will be nearer to 1° if the cosmological bound 
on neutrino masses improves further. 

The general degenerate mass matrix M in the weak basis and Qo contain too many pa- 
rameters. However, on the basis of some symmetry principle one can reduce the number 
of these parameters. One may consider a simple model of neutrino mass [8J using a lep- 
tonic Higgs doublet and three right-handed singlet fermions at TeV energy scale or below. 
Considering discrete symmetry one may obtain the degenerate mass matrix of the form 



M = m 1 (22) 




as discussed in [8]. Such pattern correspond to = and ^23 = 37r/4 in the degenerate 
mass matrix M. Then our perturbation matrix Qq in the flavor basis is 



Qo 



( e (-m + 6)/V2 {ia + h)/ \pl 

{-ia + h)/V2 {-2z + f-g)/2 {f + g)/2 | (23) 
V {ta + b)/V2 if + g)/2 {2t + f-g)/2 



This can be obtained by considering soft symmetry-breaking terms of the form mijNmNjji 
in the Lagrangian where are right-handed neutrino fields. Thus parameters mij of the 
theoretical model of neutrinos will be related to the parameters of Qq and as such Q. From 
equation. [20] it is evident that 6'i3 is very close to ~ \ rx^- ^ ^'^^ using (15) one 



gets the relative estimate of the parameters e and /. So / ~ 2 x 10^ / (me) and e is much 
smaller than that as required by CHOOZ constraint. 

Degenerate neutrino masses from an abelian family symmetry has also been considered 
in The neutrino mass matrix as considered in eq. (38) of that paper is similar in form as 
shown above in (22) & (23) provided that we consider e = and / = —g in our general form 
of the perturbation matrix Q so as to reproduce the Qo appropriate for their mass matrix. 

So far we have considered the case of different intrinsic CP parities of Majorana neutrinos 
for which there is mixing even in the exact degenerate mass limit. However there may be a 
case for which intrinsic CP parities of all the Majorana neutrinos are the same [5j. In that 
case the U matrix can be rotated away and the mixing will be controlled by the perturbation 
matrix only. So it is not possible to perform the same analysis for the perturbed part in 
the quasi-degenerate mass matrix as done in this paper. Neutrino-less double beta decay in 
such a case may yield a direct constraint on the approximate degenerate mass m, although 
^13 may not receive any additional constraint. However, same CP parity of all neutrinos is 
disfavoured from the viewpoint of the requirement of stability at the weak scale of quasi- 
degenerate neutrino masses and mixing pattern which emerge at the high scale [6] as the 
CHOOZ constraint on ^13 cannot be satisfied. In ref [H] although the same CP parity of 
quasi-degenerate neutrinos has been considered at the see-saw scale but at the weak scale it 
has been found that the appropriate mass-squared difference for oscillation of solar neutrinos 
cannot be obtained. 
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